Star formation near the Sun is driven by expansion of the Local Bubble 
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For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high- 
temperature plasma surrounded by a shell of cold, neutral gas and dust.'!* However, the precise shape and 
extent of this shell*°, the impetus and timescale for its formation®’, and its relationship to nearby star 
formation! have remained uncertain, largely due to low-resolution models of the local interstellar medium. 
Leveraging new spatial?! and dynamical constraints from the Gaia space mission”, here we report an 
analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We 
find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and 
that their young stars show outward expansion mainly perpendicular to the bubble’s surface. Tracebacks of 
these young stars’ motions support a scenario where the origin of the Local Bubble was a burst of stellar birth 
and then death (supernovae) taking place near the bubble’s center beginning ~14 Myr ago. The expansion of 
the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell 
that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing 
robust observational support for the theory of supernova-driven star formation. 


In Figure 1a (interactive), we present a 3D map of the solar neighborhood, including a new Gaia-era 3D model 
of the Local Bubble’s inner surface of neutral gas and dust!’ and the 3D shapes and positions of local molecular 
clouds constrained at ~ 1 pc resolution’. The Local Bubble’s shell is shown as a closed surface, but evidence 
suggests it could be a “Galactic chimney”, having blown out of the disk a few hundred parsecs above and below 
the Galactic plane.'* The distribution of dense gas in star-forming molecular clouds is shown with a set of 


topological “spines” derived by “skeletonizing” the clouds in 3D volume density space." 


Remarkably, we find that every well-known molecular cloud within ~200 pc of the Sun lies on the surface of the 
Local Bubble. These “surface” clouds include not just every star-forming region in the Scorpius-Centaurus (Sco- 
Cen) association (Ophiuchus, Lupus, Pipe, Chamaeleon, and Musca), but also the Corona Australis region and 
the Taurus Molecular Cloud, the latter of which lies 300 pc away from Sco-Cen on the opposite side of the bubble. 
The one exception is the Perseus Molecular Cloud, at a distance of 300 pc, which has likely been displaced by 
the recently discovered Per-Tau Superbubble"’, containing Taurus on its near side and Perseus on its far side 
(see green sphere in Figure 1b, interactive). The Taurus molecular cloud complex lies at the intersection of the 
Per-Tau Bubble and Local Bubble, displaying a sheet-like morphology consistent with being shaped by a bubble- 
bubble collision. Every Local Bubble surface cloud shows evidence of a similar sheet-like (e.g. Taurus) or 


filamentary (e.g. Corona Australis) morphology, uniformly elongated along the bubble’s surface. 


In Figure 1b, in addition to the Local Bubble surface, we show models for two kiloparsec-scale Galactic features 
discovered in the Gaia era: the Radcliffe Wave!’ and the Split". The Radcliffe Wave is a 2.7-kpc-long filament of 
gas corresponding to the densest part of the Local Arm of the Milky Way. It has the shape of a damped sinusoid, 
extending above and below the plane of the Galaxy.'* The Split, an at least 1-kpc-long gaseous feature situated 
in the disk, is argued to be a spur-like feature, bridging the Local and Sagittarius-Carina arms. Also shown 
(interactive version only) is a model of the Gould’s Belt, a disk of young stars, gas, and dust, tilted by about 20° 
with respect to the Galactic plane, which has long shaped our understanding of the architecture of the local 
interstellar medium. Prior work has suggested that the Gould’s “Belt” is a superposition of unassociated 
structures seen in projection, with all well-known regions of the "Belt" being part of either the Radcliffe Wave or 
the Split large-scale gaseous structures.'® As seen in Figure 1b (interactive), this argument about the nature of 
the Gould’s Belt is confirmed here. The right side of the assumed Gould’s Belt (the Sco-Cen association) consists 
of the entire rightward wall of the Local Bubble, while its left side consists of clouds in the Radcliffe Wave, well 
beyond the leftward wall of the Local Bubble. The Local Bubble lies at the closest distance between the Radcliffe 


Wave and the Split, with most of the dense gas at its surface co-spatial with these two kpc-scale features. 


We use measurements of the 3D positions and motions of stellar clusters to reconstruct the star formation history 
near the Local Bubble. We rely on curated samples of young stars from the literature, as summarized in Extended 
Data Table 1. Our sample includes: clusters associated with star-forming regions on the surface of the bubble 
(Taurus, Ophiuchus, Lupus, Chamaeleon, and Corona Australis), older members of the Sco-Cen association 
(Upper Centaurus Lupus, Lower Centaurus Crux, and Upper Scorpius) up to a maximum age of 20 Myr; and 
clusters in known star-forming regions along the Radcliffe Wave and the Split but beyond the boundaries of the 
Local Bubble itself (Perseus, Serpens, Orion). 


As described in the Methods section, we derive the “tracebacks” of stellar clusters associated with the Local 
Bubble and related structures. The current 3D space motions of the young stellar clusters are shown as cones in 
the interactive version of Figure 1, with the apex of the cone pointing in the direction of motion. Prior research 
has shown that the 3D space motions of the youngest clusters (S 3 Myr) can be considered probes of the 3D 
space motions of the parental gas clouds in which they were born.” Using the young stars’ motions to trace 


cloud motion, we see that not only do all star-forming clouds presently observed within 200 pc lie on the surface 


of the Local Bubble, but they also show strong evidence of outward expansion, primarily perpendicular to the 
Bubble’s surface. Tracebacks of the clusters’ motions over the past 20 Myr point to the likely origin of the Local 
Bubble — presumably the region where the supernovae driving the bubble went off. The clear implication of the 
observed geometry and motions is that all the well-known star-forming regions within 200 pc of the Sun formed 
as gas has been swept up by the Local Bubble’s expansion. 


The interactive version of Figure 1 also includes a model for Gould’s Belt!8, which illustrates that much of the 
motion previously attributed to the expansion of the assumed Gould’s Belt” is instead likely due to the 
expansion of the Local Bubble. Recent work using complementary catalogs of young stars bolster this 
interpretation, finding evidence that the Sco-Cen stellar association — a key anchor of the Gould Belt — has a 
arc-like morphology consistent with recent sequential star formation, which we now know to be triggered by 
the Local Bubble.” 


A full animation of the stellar tracebacks is provided in Figure 2 (interactive). In the static version, we show 
select snapshots at -16 Myr, -15 Myr, -14 Myr, -10 Myr, -6 Myr, -2 Myr, and the present. We observe multiple 
epochs of star formation, with each generation of stars consistent with being formed at the edge of the Local 
Bubble’s expanding shell. We find that 15-16 Myr ago, the Upper Centaurus Lupus (UCL) and Lower Centaurus 
Crux (LCC) clusters in the Sco-Cen stellar association were born about 15 pc apart, and that the Bubble itself was 
likely created by supernovae whose surviving members belong to these clusters. 


Based on the amount of momentum injection required by supernovae to sweep up the total mass of the shell 
(1.47523 x 10 ê Mo) given its present-day expansion velocity (6.7 *?3 km/s), we estimate that 15 *}’ supernovae 
were required to form the Local Bubble (see Methods section). Through an analysis of their existing stellar 
membership and an adopted Initial Mass Function (IMF), previous studies agree that UCL and LCC have 
produced 14-20 supernovae over their lifetimes.°”?! However, previous work®’ also claims that UCL and LCC 
formed outside the present-day boundary of the Local Bubble, only entering its interior in the past few 
megayears, inconsistent with an argument that they are the progenitor population. By adopting new Gaia EDR3 
estimates of the clusters’ 3D velocities, better orbit integration, and a more accurate value for the Sun’s peculiar 
motion, we find that UCL and LCC indeed coincide with the center of the bubble at its birth, lying just interior 
to its inner surface in the present day, thereby resolving this discrepancy. We explain the inconsistency between 
the stellar tracebacks for UCL and LCC proposed in this work and those from prior work in more detail in the 
Methods section.” 


Under the assumption that each star-forming molecular cloud formed due to the shell’s expansion — powered 
by UCL and LCC near its center — we fit for the temporal and radial evolution of the Local Bubble by building 
on recent analytic frameworks.” As described in the Methods section, our idealized, spherical shell expansion 
model fits for the age of the Local Bubble, the duration between supernova explosions powering its expansion, 
and the ambient density of the interstellar medium prior to the first explosion. We find that an age of 14.4*}4 


Myr, a time between supernova explosions of 1.1498 Myr, and an ambient density of 2.7 me cm” provides the 


best-fit to the dynamical tracebacks. This best-fit model for the Local Bubble’s expansion is also shown in the 


static and interactive versions of Figure 2. 


Following the presumed birth of the Local Bubble 14 Myr ago, we observe four subsequent epochs of star 
formation at the surface of its expanding shell, taking place ~ 10 Myr ago, = 6 Myr ago, ~ 2 Myr ago, and the 
present day. Around 10 Myr ago, we observe the formation of the Upper Scorpius association, as well as an 
older, recently discovered companion stellar population in Ophiuchus. Next, 6 Myr ago, both Corona Australis 
and the older stellar population of Taurus were born. Around 2 Myr ago, we detect the birth of stars in Lupus 
and Chamaeleon, as well as the younger stellar populations of Taurus and Ophiuchus. Finally, in the present 
day, we observe the current distribution of dense star-forming molecular gas, enveloping the Local Bubble. In 
Figure 2, we also overlay the solar orbit, which indicates that the Sun only entered the bubble around 5 Myr ago, 
and that it was about 300 pc away when the first supernovae went off in UCL and LCC. If this expanding shell 
scenario is true, we would expect a total of 1.74925 x 10° Mo of gas to be swept up by the Local Bubble over its 
lifetime, given our inferred ambient density of 2.7 wee cm“ and the current radius of 165 + 6 pc (see Methods 
section and Extended Data Figure 2). Based on the 3D dust currently enveloping the shell’s surface, we obtain 


an actual swept-up mass of 1.449 x 10 Mo, in agreement with the model estimate. 


The circumstances that led to the birth of the progenitor populations UCL and LCC are more difficult to 
constrain. Given the close proximity of both the Radcliffe Wave and the Split to the Local Bubble, the origin of 
UCL and LCC could be related to one of these kpc-scale gaseous structures, or to a past interaction between the 
two. While current kinematic data are limited, the 3D tracebacks of young stars in two constituent clouds along 
the Radcliffe Wave (Orion) and the Split (Serpens), but beyond the Local Bubble’s influence suggest that the 
Radcliffe Wave and the Split could have converged 20 Myr ago at the location where UCL and LCC were born. 
However, future follow-up work on the 3D motions of these two linear features would be needed to shed light 
on the true architecture of interstellar gas on kiloparsec scales at the time of UCL’s and LCC’s formation. 


Regardless of UCL and LCC’s potential origins, we find 6D (3D position and 3D velocity) observational support 
for the theory of supernova-driven star formation in the interstellar medium% — a long-invoked theoretical 
pathway for molecular cloud formation seen in numerical simulations.” The abundance of new stellar radial 
velocity data expected in Gaia DR3 should not only allow for more refined estimates of the Local Bubble’s 
evolution, but also enable similar studies farther afield, providing further observational constraints on 


supernova-driven star formation across our Galactic neighborhood. 
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Figure 1 (Interactive): A 3D spatial view of the solar neighborhood. For the best experience, please view the 
online 3D interactive version available . Panel a: A top-down projection of star-forming regions on the 
surface of the Local Bubble, whose young stars show motion mainly perpendicular to its surface. The surface of 
the Local Bubble! is shown in purple. The short squiggly colored lines (a.k.a. “skeletons”) demarcate the 3D 
spatial morphology of dense gas in prominent nearby molecular clouds". The 3D arrows indicate the positions 
of young stellar clusters, with the apex of the arrow’s cone pointing in the direction of stellar motion. Clusters 
are color-coded by longitude, as in Extended Data Table 1. The Sun is marked with a yellow cross. The zoom-in 
to the lower right shows a close-up of Ophiuchus, Pipe, Lupus, and Corona Australis on the Bubble’s surface, 
along with arrows illustrating the outward motion of their young stellar clusters. Panel b: A 3D view of the 
relationship between the Local Bubble, prominent nearby star-forming regions, and Galactic structure. The Local 
Bubble and cloud skeletons are the same as in Panel a. We also overlay the morphology of the 3D dust (gray 
blobby shapes’) and the models for two Galactic scale features — the Radcliffe Wave (red)!° and the Split (blue)"°. 
The Per-Tau Superbubble’ (green sphere) is also overlaid. The interactive version offers views from any 
direction (not just top-down), provides floating labels for star-forming regions, and includes additional layers 
(some not shown in this snapshot) which can be toggled on/off. 
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Figure 2 (Interactive): The evolution of the Local Bubble and sequential star formation at the surface of its 
expanding shell. Selected time snapshots (seen from a top-down projection) are shown here. For a full time- 
sequence, viewable from any angle (not just top-down), see the online 3D interactive version available here. The 
central figure shows the present day. Stellar cluster tracebacks are shown with the colored paths. Prior to the 
cluster birth, the tracebacks are shown as unfilled circles meant to guide the eye, since our modeling is insensitive 
to the dynamics of the gas before its conversion into stars. After the cluster birth, the tracebacks are shown with 
filled circles and terminate in a large dot, which marks the cluster’s current position. For time snapshots < 14 
Myr ago, we overlay a model for the evolution of the Local Bubble (purple sphere), as derived in the Methods 
section. The position of the local standard of rest (LSR) corresponds to the center of each panel. A more detailed 
overview of this evolutionary sequence, including the birth positions of all clusters, is provided in Extended 
Data Table 2. The solar orbit is shown in yellow and indicates that the Sun entered the Local Bubble 
approximately 5 Myr ago. 


Methods 


Deriving Stellar Cluster Properties 


In Extended Data Table 1, we summarize the properties of young stellar clusters utilized in this work. Our 
sample includes young clusters out to ~ 300-400 pc, up to a maximum age of 20 Myr. Our sample is chosen to 
provide optimal coverage of all young (<5 Myr) clusters associated with star-forming regions currently lying on 
the surface of the Local Bubble, including the Lupus”*, Ophiuchus”, Chamaeleon®, Corona Australis*, and 
Taurus Molecular Clouds”. In addition to the youngest stellar populations on the Local Bubble’s surface, several 
recent studies provide evidence that at least two of the surface clouds — Taurus? and Ophiuchus” — exhibit 
multiple generations of star formation, with an older stellar population (> 5 Myr) existing alongside the younger 
one. For these clouds, we include both young and old populations. Older stellar populations in the Scorpius 
Centaurus*4 association —- Upper Centaurus Lupus, Lower Centaurus Crux, and Upper Scorpius are likewise 
included in Extended Data Table 1. While many of the known moving groups (e.g. Beta Pic, Octans, and Carina?) 
also lie inside the bubble, their ages (> 25 Myr) are larger than the bubble’s estimated age and are thus excluded 
from our analysis®. Finally, we also include clusters beyond the Local Bubble but associated with star-forming 


regions along the Radcliffe Wave and the Split, including the Perseus*, Orion”, and Serpens% Molecular Clouds. 


We rely on existing studies (see Extended Data Table 1) to determine stellar membership of each cluster, but we 
only include stars which are detected in the Gaia EDR3 catalog. We uniformly associate stellar members of each 
cluster with Gaia EDR3 using a crossmatch radius of 1” to obtain its sky coordinates, parallax, and proper 
motions. If radial velocities for the stars are provided along with the target selection in their original publication, 
we adopt those existing radial velocities in our analysis, some of which are obtained with high-quality ground 
based near-infrared spectroscopy. Otherwise, if no radial velocity data are provided, we utilize the Gaia radial 
velocities and restrict our analysis to only those stars which also have a Gaia radial velocity detection. We largely 
rely on the sample selection outlined in each cluster’s original publication to filter outliers, many of which are 
defined using Gaia DR2 data. However, as an additional constraint, we require that all stars must have a parallax 
over error greater than two and small renormalized unit weight error (RUWE < 1.4).*” We also require the radial 
velocity error to be < 5 km/s, a relatively generous cut chosen because our algorithm incorporates the errors on 
the individual stellar measurements when determining the mean cluster motion. After applying all these cuts, 
we perform a sigma-clipping procedure using the astropy*® package, removing extreme outliers whose radial 
velocities are inconsistent with the rest of the cluster population at the 3a level. Finally, we require that each 
cluster have at least three stellar members. The mean stellar membership is much higher than this, averaging 37 


stars per cluster. 


In order to transform the sky coordinates (right ascension a and declination ò), parallax (7), proper motions 
(ura, unec), and heliocentric radial velocities (Vnelic) of members to an average 3D space position and 3D velocity 


of the cluster in a Heliocentric Galactic Cartesian reference frame (x, y, z, U, V, W), we utilize the extreme 


deconvolution algorithm.” The Extreme Deconvolution algorithm infers an n-dimensional distribution function 
using a Gaussian Mixture Model (GMM) given the presence of noisy, incomplete, and heterogeneous samples 
of the underlying population. We apply the algorithm to infer the average 6D phase information of each cluster, 
given the observed values and estimated error covariances of its stellar members. For each star we use the 
astropy** functionality within galpy*® to compute the star’s Heliocentric Galactic Cartesian coordinates (x, y, z) 
and associated space motions (U, V, W) given the observed Gaia quantities (a, 5, Tl, URA, LDEC, Vrelio). We assume U 
points towards the Galactic center, V points toward the direction of Galactic rotation, and W points toward the 
North Galactic Pole. In order to accurately estimate errors on (x, y, z, U, V, W) for each star, we randomly sample 
one hundred times from a multi-dimensional Gaussian in (a, 6, T, URA, MDEC. Vhelio) Space assuming Gaussian 
uncertainties on all parameters as reported in the Gaia EDR3 catalog. Transforming each sample to a Heliocentric 
Galactic Cartesian coordinate frame, we then calculate the covariance of the set of samples for each star. Feeding 
the set of (x, y, z, U, V, W) values for the individual stellar cluster members and their associated sample 
covariances into the Extreme Deconvolution algorithm, we obtain the mean and variance of a single 6D Gaussian 
defining the average 3D position and 3D space motion of each cluster, as shown in Extended Data Table 1. We 
adopt a peculiar solar motion of (Uo, Vo, Wo) = (10.0, 15.4, 7.8) km/s*! and correct the (U, V, W) values of each 
cluster for this solar motion to obtain its current 3D space velocity with respect to the Local Standard of Rest 
(LSR) frame (Utsr, Visr, Wiser). We use the mean cluster motion to “traceback” its trajectory in the Galaxy. To 
“traceback” a cluster means to compute the 3D position and 3D motion that the cluster would have had at 
different times in the past, given estimates of its present 3D position and 3D motion. In practice, the full bound 
orbit of the cluster can be computed, and the small section of the orbit constituting its trajectory in the recent 


past is extracted. 


We perform the dynamical traceback of each cluster using the galpy package”, which supports orbit integrations 
in a Milky Way-like potential, consisting of a bulge, disk, and dark matter halo. We sample the orbit from -20 
Myr to the present day. Alongside each cluster, we also trace the Sun’s orbit backward in time over the past 20 
Myr and correct each cluster’s orbit for the Sun’s peculiar motion to ensure all orbits remain in the Local 
Standard of Rest frame. We emphasize that galpy only accounts for the gravitational potential of the Galaxy and 
does not consider the gravitational effects of the parental gas clouds in which many of these clusters are 
embedded. However, given the large extent of the Local Bubble relative to any individual star-forming region, 
and the fact that the Galactic potential should dominate over the gravity of any local gas, galpy still serves as a 
useful probe of the dynamics of the Local Bubble. This argument is bolstered by recent results from numerical 
simulations, which indicate that stellar orbits can be recovered with high fidelity up to 20 Myr in the past, even 
without explicitly modeling non-axisymmetric components of the potential.“ The dynamical tracebacks for all 
clusters in Extended Data Table 1 are publicly available at the Harvard  Dataverse 
(https://doi.org/10.7910/DVN/E8PQOD). 


Modeling the Local Bubble’s Expansion 


In this section, we derive an analytic model for the radius and expansion velocity of the Local Bubble as a 


function of time, using constraints provided by the dynamical traceback data summarized in Extended Data 


Tables 1 and 2 (see the Data Availability section to access the full traceback data on each cluster). The results of 
this section underpin our model for the temporal evolution of the Local Bubble shown in the static and 


interactive versions of Figure 2. 


To model the expansion of the Local Bubble, we utilize recent literature that leverages 1D spherically-symmetric 
hydrodynamic simulations using the Athena++ code to study the dynamical evolution of superbubbles driven 
by clustered supernovae in a uniform medium.” Specific treatment is given to the effects of cooling at the 
shell/bubble interface. Building on this recent literature, we adopt an analytic model” describing the radius, R, 
of the superbubble’s expanding shell as a function of time t since its birth, parameterized by the ambient density 
no of the interstellar medium, the cooling efficiency at the shell’s surface 0, the time separation between 
supernovae explosions Atsne within the cluster powering its formation,, and the energy input per supernova 


explosion Esn, as follows: 


=A. -1 3 


R(t) = 83 pe x (1-0)! (F (Ate) (2) (1) 0) 


1051 erg 0.1 Myr 1 cm7 1 Myr 


At any time t the (x, y, z) coordinates of the surface of this expanding spherical shell, centered on (Xcen, Ycen, Zcen), 


corresponding to the epicenter of the supernova explosions, can be parameterized as: 


(x(t) = Niels + (y(t) ~ View) + (z(t) = Zeen)” = R(t)? (2) 


The 3D positions of each cluster at birth (derived from the dynamical tracebacks given the cluster’s age) 
provide a constraint on the bubble’s radius as a function of time. So, under the assumption that the formation 
of the young stars listed in Extended Data Table 2 was triggered by the shell’s expansion, we can infer the 
parameters governing the Local Bubble’s evolution using this analytic framework. We emphasize that this 
theoretical formalism is an approximation of the bubble’s true morphology. We do not actually expect the 
bubble to expand spherically because the interstellar medium is highly turbulent with significant density 
fluctuations. Indeed the Local Bubble today is observed to have a complex, non-spherical morphology. 
Extended Data Table 2 lists the 3D (x, y, z) birth positions for the subset of young clusters utilized to model the 
bubble expansion in the Local Standard of Rest Frame, given their stellar tracebacks and estimated ages. 


Since several of the parameters governing the superbubble’s evolution are degenerate, we make a number of 
simplifying assumptions. Prior work®7?! estimates that 14-20 supernovae in Upper Centaurus Lupus (UCL) 
and Lower Centaurus Crux (LCC) stellar groups over the past ~13 Myr have together created the Local Bubble. 
The formalism of the analytic superbubble model assumes that all supernovae are driven from a single 
location. Rather than model each individual supernova explosion given the tracebacks of UCL and LCC, we 
assume that the epicenter of the explosion (Xcen, Ycen, Zcen) lay equidistant between UCL and LCC at the time of 
the first explosion, texp. This approximation is justified as UCL and LCC lay roughly co-spatial at early times, 
when the most powerful supernovae driving the superbubble’s expansion would have gone off. We leave the 


time of the first supernova explosion, tep, as a free parameter in our model. The subsequent evolution of the 


Local Bubble is governed by tex and Equation (1). We assume a fixed energy input per supernova of 10°! erg. 
We also assume a fixed cooling efficiency 0 of 0.7. However, we test a variety of cooling efficiencies, ranging 
from 0.4 - 0.9, and find that fixing the cooling efficiency to a value of 0.7 does not affect our estimate for the 
time of the first explosion, and only has a modest effect on the ambient density and duration between 
supernovae (with the variation falling within our reported uncertainties on these parameters, as we will later 


show in Extended Data Figure 1). 


Having fixed the cooling efficiency and energy input per supernova , the free parameters of our model include 
the ambient density, no, the time between supernova explosions, AtsNe, and the time of the first supernova 

explosion tep. We infer the values of no, Atsne, and tepin a Bayesian framework. We assume that the density of 
the shell has a Gaussian profile with an uncertainty (or thickness) of Ar, which corresponds to a log-likelihood 


of the following form: 


Log(L) = = ry, [in(2n Ar VE [R(ti no, Atayo) - tori (3) 
Here, the R(ti, no, Atsne) term is the radius of the Local Bubble’s expanding shell governed by Equations (1) 
and (2) , evaluated at time ti, corresponding to the difference between the time at which the ith cluster was 
born and the time of the first explosion (ti= tbirni- tep). The Local Bubble’s shell with radius R(ti) is centered on 
(Xcen, Ycen, Zcen), derived from the mean 3D position of UCL and LCC in the Local Standard of Rest at time tep, 
when the first supernova went off in UCL or LCC. The r(ti) term is the radius of a sphere with the same center 
as R(ti), such that the ith cluster born at traceback time ti lies on its surface, with coordinates of: 


(x(t;) = Xien)” + (y (ti) = Veen)” + (z(t) = Zn = r(t)? (4) 


Finally, the Ar term in the log-likelihood given in Equation (3) should be interpreted as an error term, 
encompassing uncertainties in both the ages of the clusters and on their mean (Ursr, Visr, Wisk) motions. We 
infer Ar as an additional free parameter in our model. The total log-likelihood of all n clusters is the sum of 
their individual log-likelihoods, evaluated at the respective time of their births, which will be optimized when 
the difference between R(ti) and r(ti) is minimized. 


Using the log-likelihood in Equation (3), we sample for the values of fexp, no, AtsNe, and Ar using the nested 
sampling code dynesty.* For tep, we adopt a truncated normal prior with a mean of -13 Myr and a standard 
deviation of 1 Myr over the range -16 Myr to -10 Myr (consistent with previous evolutionary synthesis models 
of UCL and LCC).°?! For no, we adopt a truncated log-normal prior with a mean of 2 cmand a standard 
deviation of a factor of two over the range 0.1 to 10 cm*, consistent with the density range explored in the 
Athena++ simulations underpinning the expansion model.” For AtsNe, we adopt a truncated log-normal prior 
with a mean of 0.8 Myr and a standard deviation of a factor of two, over the range 0.05 Myr to 3 Myr 
(consistent with previous estimates of approximately 16 supernovae occurring in UCL and LCC over the past 


13 Myr)’. Finally, guided by the typical errors on the 3D motions (a few km/s) and the ages (a few Myr), we 


adopt a truncated normal prior on Ar, with a mean of 15 pc and a standard deviation of 5 pc, over the range 0 


to 30 pc. We run with the default parameters of dynesty’s dynamic nested sampler. 


The results of our sampling procedure are summarized in Extended Data Figure 1. We find a median value 
and 1o errors (computed using the 16th, 50th, and 84th percentiles of the samples) of Atsne= 1.0619$3 Myr , no= 
2.7 TE cm, texp = —14.39+0:78 Myr, and Ar= 23.31133 pc. The best-fit model corresponds to an epicenter of the 
bubble of (xcen, Ycen, Zeen) = (39, 7, -18) pc in the LSR frame 14.4 Myr ago. Adopting these median parameters, a 
model for the evolution of the Local Bubble is overlaid in Figure 2. We compute the radius of the Local 
Bubble’s shell and its expansion velocity over its lifetime, as plotted in Extended Data Figure 2 (panels a. and 
b.). Based on our model, and leveraging the full set of posterior samples, we estimate a present-day expansion 
speed of 6.7 +93 km/s and a radius of 165 + 6 pc. This present-day expansion speed is consistent with the 
current range of 3D velocity magnitudes of stars at the surface of its shell (~ 5 - 9 km/s), assuming that the rest 
velocity of the shell lies within a few km/s of the Local Standard of Rest (as we expect it to, since the LSR is 


currently inside the Bubble). 


However, as seen in Extended Data Figure 1, there is a strong covariance between Atsne and nosuch that an 
increase in the density of the ambient medium can be compensated for by a decrease in the time between 
supernova explosions, and vice versa. Given the limitations of our modeling, we intend the superbubble 
evolution shown in Figure 2 to only serve as a possible, idealized, example of how the Local Bubble could have 


reached its present-day morphology. 
Potential Origin of the Local Bubble 


In this section, we seek to shed additional light on the origin of the Local Bubble, by using a momentum analysis 
to test whether UCL and LCC harbored enough supernovae to excavate a cavity the size of the Local Bubble. To 
do so, we compute the momentum of the Local Bubble’s shell from its current expansion velocity (calculated in 
the previous section and shown in Extended Data Figure 2) and estimates of its mass (obtained from 3D dust 
mapping). We can then further build on the analytic superbubble model constrained above to obtain the 
expected average momentum injection per supernova, Ê. The ratio of the total momentum of the shell to the 
average momentum injected per supernova provides a constraint on the number of supernovae required to 
power its expansion, which can then be compared with existing estimates for how many supernovae have gone 


off in UCL and LCC based on population synthesis modeling.”7! 


To obtain the momentum of the shell, we calculate its swept-up mass Msheu by integrating the 3D volume density 
derived from the 3D dust map’ between a distance of [Rsnen, Rshen + Rtnickness] from the Sun and multiplying by the 
mean particle mass (m) = 1.4m, where muis the mass of a proton, and the factor of 1.4 corrects for the helium 
abundance. Rsnet is the boundary (i.e. inner radius) of the Local Bubble (shown in Figure 1) and Rinickness is the 
shell’s thickness, such that Rshet' + Rtnickness corresponds to the outer radius. The Local Bubble model we utilize 
estimates Rtnickness to be between 50 - 150 pc!, which is quite large, but encompasses the full depth of structure 


currently lying on the bubble’ surface. 


Adopting Rinickness = 100 pe with an estimated 1ø uncertainty of 50 pc, we obtain Msnei = 1.44983 x 10° Mo for the 


swept-up mass. To estimate the current expansion velocity of the Bubble vexp, we leverage the posterior samples 
from our model describing the Local Bubble’s evolution (fit to the dynamical tracebacks; see Extended Data 
Figure 1) to evaluate the velocity of the shell vsnen = dR/dt in the present day (see Extended Data Figure 2). Doing 
so, we obtain vep = 6.7 +?3 km/s. To propagate uncertainties throughout this section, we use the full set of 
posterior samples when leveraging parameters from the expansion model. To propagate uncertainties on the 
observed swept up mass of the shell, we draw the same number of samples from a Gaussian with a mean of 
Mstett= 1.4 x 10° Mo and a standard deviation of 6 x 105 Mo. Whenever uncertainties are reported, we use the 
16th, 50th, and 84th percentiles of the samples to compute the median and 1o error bounds. Given samples for 


the current expansion velocity and swept-up mass, the corresponding momentum of the Local Bubble’s shell is: 
Pshell = Mosneni x Vexp = 9.6 ar x 10 6 Mo km/s 


In the previous section, we utilized an analytic superbubble expansion model” to constrain the temporal 
evolution of the Local Bubble’s size and velocity. The same simulations underpinning this formalism indicate 
that the average momentum injected per supernovae 6 depends on the same free parameters, namely the energy 
input per supernova explosion, Esn, the cooling efficiency, 8, the duration between supernova explosions, Atsne, 


the density of the ambient medium, no, and time t: 


1 


p(t) = 4 x 10° Mo kms71 (1 — 9)s (i) (2e) (- a ( Ly 5 


Therefore, we again fix 0 = 0.7, Esn = 10°! erg, and leverage the samples of Atsne, no, and terp. We evaluate 6 over 


the lifetime of each realization of the bubble and draw a random sample of 6 from each distribution. We then 

take the 16th, 50th, and 84th percentiles of the random samples drawn from all realizations (see Extended Data 

Figure 1) to obtain a mean value and 1ø uncertainties on the average momentum injected per supernova £: 
p=65 +46 x 105 Mokm/s 

Finally, dividing all samples of psnen by all samples of 6 we obtain: 


None = psnet / = 15 +4’ supernovae 


The marginal distribution of Nsne is shown in Extended Data Figure 3. The average number of 15 +1 supernovae 
is in good agreement with previous results®””!, which argue that UCL and LCC have produced 14-20 supernovae 


over their lifetimes based on an analysis of their current stellar membership and modeling of a Salpeter“ IMF. 


Consistent with this physical scenario, before supernovae start going off in UCL and LCL, our model predicts a 
typical ambient density of the interstellar medium of no= 2.7 A cm. Assuming that the volume of sphere 
with a radius of 165 + 6 pc (our current estimate for the radius of the Local Bubble from the analytic expansion 
model) is uniformly filled with a gas density of no= 2.717733 cm*, we would expect 1.7725 x 10 ê Mo of gas to 
be displaced and swept up into its surrounding shell over its lifetime. Using the 3D dust measurements, we 
measure an actual swept-up mass of Msnen = 1.4063 x 10° Mo on the surface of the Local Bubble, aligned with 


this estimate of 1.7 million Mo. 


Given uncertainties in the exact volume density of gas and dust inside the bubble traced by 3D dust - with some 
recent work suggesting that there is only a modest drop in volume density coincident with the cavity of hot 
ionized gas — we note that these momentum calculations (and more broadly the analytic superbubble 
expansion model) are largely insensitive to the exact difference in density interior and exterior to the bubble’s 
boundary. The analytic expansion model we adopt” includes the ambient density of the interstellar medium at 
the time of the first supernova explosion as a free parameter; however, this ambient density parameter is 
constrained using the dynamical traceback data (not the 3D dust maps) and does not require any explicit 
assumptions about the exact density interior and exterior to the bubble as the shell expands. Similarly, while the 
swept-up mass calculation does depend on 3D dust mapping’, the mass is only measured within the denser shell 
of neutral gas and dust (where the uncertainties on the underlying dust extinction are smaller) and has no 
dependence on the density in the hot inner cavity, which could be subject to larger uncertainties due to the very 


low levels of dust extinction. 
The Stellar Tracebacks of UCL/LCC and the Peculiar Motion of the Sun 


Tracebacks of the average stellar cluster motions are expressed relative to the “Local Standard of Rest,” which 
itself depends on the measured motion of the Sun. Even in the age of Gaia, the Sun’s exact motion in the Galaxy 
is uncertain, especially along the “Y” direction (along Galactic rotation). In this paper, we revise estimates of 
the tracebacks of the UCL and LCC clusters from prior literature®’, which argues that UCL and LCC were born 
outside of the present-day boundary of the Local Bubble even though their members were the likely progenitors 
of the Local Bubble. In this section, we revisit these extant stellar traceback calculations from prior work®’ finding 


that a revised calculation places UCL and LCC near the center of the Local Bubble when they formed. 


While previous studies”! have proposed UCL and LCC as the likely progenitor population, the prior literature 
in question’ is the first study to perform an unbiased stellar search of all nearby B stars to track down the remains 
of OB associations hosting supernovae capable of powering the Local Bubble. This prior literature® confirms that 
UCL and LCC are the only populations capable of having powered the Local Bubble, but, after tracing back the 
stellar members, find that both UCL and LCC formed > 100 pc outside the present-day boundary of the Local 
Bubble. Extended Data Figure 4a shows the extant stellar traceback results, where UCL and LCC only lay interior 
to the present-day boundary of Local Bubble’s during the past few megayears.° The prior literature® noted this 
potential inconsistency (how does a supernova cause a bubble it does not lie within?), but argue that the location 
of the UCL and LCC clusters with respect to the center of the bubble is not crucial. 


The prior literature’ which places UCL and LCC outside the bubble use Hipparcos parallaxes and proper 
motions with radial velocities collected from the literature to derive the (U, V ,W) 3D space motion of each 
member. Using the (U, V ,W) value of each individual star, prior literature® calculates the mass—weighted mean 
(U, V,W) velocity of all members of UCL and LCC, and then assign this group 3D U, V, W space velocity to each 
individual stellar member for the purposes of performing the tracebacks. The prior literature® also states that 
they correct the mean 3D group velocity for the Sun’s peculiar motion, adopting a value of (Uo, Vo, We) = (10.0, 
5.2, 7.2) km/s“, so that the traceback of each star is reported in the Local Standard of Rest frame. 


As a first step, we attempt to reproduce the original results from prior literature’ (see Extended Data Figure 4a) 
using their own data. Specifically, we calculate the mass-weighted mean velocity of the UCL and LCC stars from 
their Table A16 and Equation 2°, finding (U, V, W = -7.1 -20.6, -5.8) km/s without correcting for the solar motion. 
With a value for the Sun’s peculiar motion of (10.0, 5.2, 7.2) km/s* used in the prior literature, these values 
translate to (Utsr, Visr, Wisr = 2.9, -15.4, 1.4) km/s in the LSR frame. Extended Data Figure 4b shows dynamical 
tracebacks we calculate in galpy (see Methods) from the Table A1 data in prior literature and with the value for 
the Sun’s peculiar motion of (10.0, 5.2, 7.2) km/s*. Using the default orbit integrator in galpy with its standard 
Milky Way potential”, we are unable to reproduce the results from prior literature, particularly the strong curvature 
in the tracebacks along the +X direction. The prior literature uses an epicyclic approximation for the stars’ 
motions (see their Section 2)°, which may be responsible for part of the discrepancy. This systematic smearing 
out of the tracebacks toward +X in their Figure 2° manifests in the entire sample of B stars and is not isolated to 
UCL and LCC. 


In examining the stellar motions of UCL and LCC, we find that the choice of the solar peculiar motion — necessary 
to convert to the LSR frame — has a non-negligible effect on where the birthplaces of UCL and LCC lie with 
respect to the center of the Local Bubble. There have been dozens of attempts to constrain the Sun’s peculiar 
motion, but Vo (motion in the direction of Galactic rotation, along Y) remains highly uncertain: current estimates 
range between Vo = 4 and 16 km/s.” 


The value of Vo = 5.2 km/s** adopted in prior literature®’ is one of the lowest values of Vo. As Extended Data 
Figure 4c shows, if we use the same Hipparcos data utilized in the prior literature but replace their value for Vo 
(5.2 km/s**) with a value of 15.4 km/s*! we find that UCL and LCC are born near the bubble’s center. Extended Data 
Figure 4d shows that updating the Hipparcos data with Gaia data and adopting the same peculiar motion used 
in this work": makes almost no difference to the tracebacks, suggesting that uncertainty in Sun’s peculiar motion 
is the dominant source of uncertainty in determining the birth location of UCL and LCC. We adopt a value of 
15.4 km/s*! when calculating all the tracebacks in this paper, toward the upper end of estimates for Vo. Given 
the large uncertainty on Vo, we have tested the robustness of our results to the choice of solar peculiar motion 
and find that any motion for Vo 2 10 km/s is entirely consistent with the physical scenario we propose here. This 
Vo=10 to 16 km/s range encompasses the vast majority of estimates for Vo used in the field today. 


Data Availability 


The datasets generated and/or analyzed during the current study are publicly available on the Harvard 
Dataverse (https://dataverse.harvard.edu/dataverse/local bubble_star_formation/), including Extended Data 
Table 1 (doi:10.7910/DVN/ZU97QD), Extended Data Table 2 (doi:10.7910/DVN/1VT8BC), per-star data for 
individual stellar cluster members (doi:10.7910/DVN/1UPMDX), and the cluster tracebacks 
(doi:10.7910/DVN/E8PQOD). 


Code Availability 


The results generated in this work are based on publicly available software packages and do not involve the 
extensive use of custom code. Given each star’s reported Gaia data, we use the astropy® package to obtain the 
Heliocentric Galactic Cartesian positions and velocities. The extreme deconvolution algorithm in the astroML** 
package is used to estimate the mean 3D positions and velocities of the stellar clusters. The Orbit functionality 
in the galpy* package is used to perform the dynamical tracebacks. The dynesty® package is used to fit the 
analytic superbubble expansion model and determine the best-fit parameters governing the Local Bubble’s 


evolution. 
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Cluster Region Cluster Subgroup Sample Nstars | b d x ox y oy Z oz Uisr usr Visk Ouse Wisr owsr Age 


si ° pce pe pe pe pe pe pe km/s km/s km/s km/s km/s km/s Myr 


(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) 
Perseus NGC1333 Ortiz-Leon etal. 2018 28 158.3 -20.5 293 -255 9 101 3 -102 3 -72 09 50 12 -21 09 1 
Perseus 1C348 Ortiz-Leon et al. 2018 126 160.5 -17.8 315 -2822 8 100 3 -96 3 69 10 96 09 01 09 3 
Taurus C2-L1495 Krolikowski etal. 2021 24 168.9 -15.8 130 -122 2 24 1 -35 2 -6.0 05 32 07 -31 05 1 
Taurus C8 - B213 Krolikowski etal. 2021 12 170.9 -15.8 159 -151 4 24 2 43 2 79 03 24 05 07 02 3 
Taurus D4 - North Krolikowski et al. 2021 21 172.3 -16.3 135 -128 25 17 7 -38 9 -49 29 52 59 -14 31 2 
Taurus C6 - L1524 Krolikowski etal. 2021 32 173.7 -15.5 128 -123 2 13 4 -34 1 -58 09 36 16 -19 09 2 
Taurus C7 - L1527 Krolikowski etal. 2021 11 174.0 -13.6 140 -136 3 14 1 -33 1 -59 13 40 07 -22 07 3 
Taurus C5 -L1546 Krolikowski et al. 2021 13 175.6 -16.5 162 -155 3 12 1 -45 1 -67 05 17 04 06 07 2 
Taurus D3 - South Krolikowski etal. 2021 13 176.7 -18.0 122 -116 3 6 10 -37 7 -414 18 76 26 -20 08 6 
Taurus C1-L1551 Krolikowski etal. 2021 14 179.1 -20.0 145 -136 2 2 1 -49 1 -54 17 03 07 05 11 2 
Taurus D2 - L1558 Krolikowski etal. 2021 13 179.9 -20.0 131 -123 29 0 8 -45 11 -50 25 43 52 -14 14 3 
Taurus D1 - L1544 Krolikowski etal. 2021 12 181.3 -7.3 172 -170 8 -3 5 -21 6 -89 06 15 147 40 12 3 
Taurus C9 - 118TauEast Krolikowski etal. 2021 4 184.0 -4.9 109 -109 2 7 1 9 2 -34 04 -34 02 -12 01 6 
Orion L1616 Großschedl etal. 2021 5 203.6 -24.6 382 -318 15 -139 6 -159 6 -88 27 47 11 16 10 2 
Orion NGC2068/2071 Großschedl etal. 2021 44 205.2 -14.3 413 -362 10 -170 5 -102 3 -13.6 13 32 10 -05 13 2 
Orion NGC 2023/2024 Großschedl etal. 2021 37 206.6 -16.3 401 -344 10 -172 5 -112 3 -135 18 26 16 05 12 2 
Orion Sigma Ori Großschedl etal. 2021 38 206.8 -17.3 401 -342 7 -173 4 -119 3 -16.3 10 -05 08 04 11 3 
Orion NGC1977 GroRschedl etal. 2021 42 208.5 -19.1 393 -326 6 -177 3 -128 2 -151 13 -0.9 10 -11 11 3 
Orion Orion A, Head Großschedl etal. 2021 241 209.2 -19.5 393 -324 9 -181 5 -131 4 -122 20 16 19 08 16 2 
Orion Orion A, Tail Großschedl etal. 2021 153 211.8 -19.4 408 -327 17 -203 19 -135 8 -8.2 27 28 15 05 10 2 

Chamaeleon 1-North Galli et al. 2021 34 297.0 -14.8 192 84 1 -165 2 -49 1 -06 06 -38 07 -39 06 2 
Chamaeleon 1-South Galli et al. 2021 34 297.2 -15.8 188 82 2 -161 3 -51 2 16 06 -42 11 -34 06 2 
Sco-Cen Towar GOEL Gns, Gagne etal. 2018a,b 42 301.3 6.0 110 57 14 -933 13 11 15 19 15 61 20 19 16 15 
Chamaeleon 2 Galli et al. 2021 7 303.7 -14.4 196 105 2 -158 2 -48 2 -08 19 -32 33 -08 14 2 


Sco-Cen la sae Gagne et al.2018a,b 40 330.7 13.3 129 109 25 -61 17 29 16 48 35 46 19 27 14 16 


Lupus 4 Galli et al. 2020 5 3364 84 160 145 1 -63 1 23 1 5.2 15 -20 09 02 03 3 
Lupus 3 Galli et al. 2020 27 3395 94 158 146 2 -54 1 25 1 67 27 -20 13 02 09 3 
Lupus Off-Cloud Population Galli et al. 2020 3 340.1 109 158 146 4 -52 1 29 5 58 36 -14 12 02 05 3 
Sco-Cen Upper Scorpius, USCO Gagne etal. 2018a,b 43 351.4 22.0 142 130 8 -19 7 53 9 51 21 -08 12 09 1.7 10 
Ophiuchus Rho Oph, Population | Grasser et al. 2021 96 353.1 17.2 138 131 3 -15 1 40 2 4.7 1.7 0.3 09 -14 1.0 1 
Ophiuchus Rho Oph, Population II Grasser et al. 2021 32 353.9 18.3 139 131 6 -14 2 43 3 §9 22 04 13 18 08 10 
Corona Australis Off-Cloud Population Galli et al. 2020 12 358.3 -13.7 146 141 3 -4 3 -34 4 39 22 -22 04 06 1.0 6 
Corona Australis On-Cloud Population Galli et al. 2020 6 359.9 -16.7 149 143 3 0 1 43 2 46 21 -144 04 -140 06 5 
Serpens Far South Herczeg et al. 2019 3 273 24 384 341 3 176 0 14 0 123 26 19 17 33 02 3 


Extended Data Table 1: Summary of the 3D positions and 3D velocities of young stellar clusters within 400 
pc of the Sun. The rows are color-coded according to the cluster’s longitude and map to the colors of the clusters 
shown in Figures 1 and 2. (1) Name of the cluster region. (2) Name of the cluster subgroup. (3) Original 
publication forming the basis for stellar membership. (4) Number of stellar members used to estimate the 
average cluster properties. (5) Mean longitude (6) Mean latitude (7) Mean distance (8-9) Mean x position and x 
uncertainty in Heliocentric Galactic Cartesian coordinates. (10-11) Mean y position and y uncertainty in 
Heliocentric Galactic Cartesian coordinates. (12-13) Mean z position and z uncertainty in Heliocentric Galactic 
Cartesian coordinates. (14-15) Mean U motion and U uncertainty along the Heliocentric Galactic Cartesian x 
direction in the LSR frame. (16-17) Mean V motion and V uncertainty along the Heliocentric Galactic Cartesian 


y direction in the LSR frame. (18-19) Mean W motion and W uncertainty along the Heliocentric Galactic 


Cartesian z direction in the LSR frame. (20) Estimated age of the cluster drawn from the publication in column 


(3). 


Timestep (Myr) 


-16 -15 -14 -10 -6 -5 -3 -2 -1 Now 
Geran USCO and CrA, Taurus “Young”, Taurus “Young”, Taurus “Young”, Dense Gas 
UCL Born LCC Born Make Oph “Old” Taurus “Old” CrA Born Lupus Cham Oph “Young” Enveloping 
Born Born Born Born Born Bubble 
Bubble 

UCL (40,19,-24) (44,14,-20)  (48,9,-17) (64,-11,-3) (81,-32,10) (86,-37,14) (95,-47,20) (99,-51,23) (104,-56,26)  (109,-61,29) 
Lec (31,6,-18)  (32,0,-16) (39,-27,-9) (45,-54,0) (47,-61,1) (51,-74,5) (53,-80,7) (54,-87,9) (57,-93,11) 

Rho Oph “Old” 
‘Pocuslion i) (75,-8,14) (97,-11,28) (102,-11,31) (113,-12,37)  (119,-13,39)  (125,-13,41)  (131,-14,43) 
— (83,-10,29) (100,-14,42) (105,-15,44) (115,-17,49)  (119,-18,50) (125,-18,52) (130,-19,53) 


Taurus 


(D3 - South) (-92,-40,-22) (-96,-32,-25) (-103,-16,-30)  (-107,-8,-33) (-111,-1,-35)  (-116,6,-37) 


(of. Cloud Population) (120,9,-34) (123,7,-34)  (130,2,-35)  (134,0,-35)  (137,-1,-35)  (141,-4,-34) 


(co- ered (-89,13,-1)  (-92,9,-2) (-98,2,-5) (-101,0,-6) (-104,-4,-8)  (-109,-7,-9) 


Corona Australis 


(On-Cloud Population) (120,7,-34)  (129,4,-38) (133,2,-40) (138,1,-41) (143,0,-43) 


— (129,-56,22)  (134,-59,22)  (139,-61,23)  (145,-63,23) 


. 
(c7-L1527) (-117,1,-25)  (-123,6,-28)  (-129,10,-30) _(-136,14,-33) 


a 
(C8 - B213) (-127,16,-44)  (-135,19,-44)  (-143,21,-43) (-151,24,-43) 


= (126,-48,24)  (132,-50,25)  (139,-52,25) — (146,-54,25) 


Taurus 
(D2 - L1558) (-108,-12,-39)  (-113,-8,-41) (-118,-4,-43)  (-123,0,-45) 


z 
(01 211544) (-143,-8,-18)  (-152,-6,-19)  (-161,-5,-20) (-170,-3,-21) 


Taurus 
(C6 - L1524) (-110,6,-29) (-116,9,-32)  (-123,13,-34) 


Erri (86,-152,-43) (84,-156,-47) (82,-161,-51) 


Taurus 
(C5 - L1546) (-141,8,-46) (-148,10,-46) (-155,12,-45) 


Taurus 
(D4 - North) (-118,6,-34) (-123,12,-36) (-128,17,-38) 


Taurus 
(C1 - L1551) (-126,1,-50) (-131,1,-49) (-136,2,-49) 


Oana (85,-157,-40)  (84,-161,-44)  (84,-165,-49) 


— (107,-151,-46)  (106,-154,-47) (105,-158,-48) 


Taurus 
(C2 - L1495) (-116,20,-32) (-122,24,-35) 


A PE SE (128,-48,28)  (134,-49,29)  (140,-51,29) (146,-52,29) 


Rho Oph “Young” 


(Population 1) (126,-16,42)  (131,-15,40) 


Extended Data Table 2: Temporal evolution of cluster births at the surface of the Local Bubble’s expanding 
shell. The time sequence from left to right shows the (x,y,z) positions (in pc) of clusters as a function of time. If 
a cell appears blank, the cluster has not yet been born, with the first filled cell (from left to right) indicating the 
birth position of the cluster. These birth positions are used to fit the analytic superbubble expansion model 
presented in the Methods section. All positions are given with respect to the Local Standard of Rest. The rows 
are color-coded according to the cluster’s longitude and map to the colors of the clouds shown in Figures 1 and 


2, as well as Extended Data Table 1. Significant events corresponding to each time step are summarized at top. 
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Extended Data Figure 1: 1D and 2D marginal distributions (“corner plot”) of the model parameters 
governing the evolution of the Local Bubble’s expanding shell. Parameters include the time since the first 
explosion (i.e. the age of the Local Bubble), tep, the density of the ambient medium the bubble is expanding 
into, no, the time between supernova explosions powering its growth, Atsne, and the thickness/uncertainty on 
the expanding shell radius Ar. In the 1D distributions, the vertical dashed lines denote the median and 10 


errors, while in the 2D distributions, we show the 0.50, 1a, 1.50, and 20 contours. 
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Extended Data Figure 2: Temporal evolution of the Local Bubble, based on the fit to the dynamical tracebacks 
and the analytic expansion model? summarized in the Methods section. Panel a) The evolution of the Local 
Bubble’s expansion velocity vep. Panel b) The evolution of the Local Bubble’s shell radius Rsnen. Panel c) The 
evolution of the average momentum injection per supernova /. The thick purple line represents the median fit, 
while the thin purple lines represent random samples. We estimate a current radius of 165 + 6 pc and current 


expansion velocity of 6.7 +89:37 km/s, corresponding to time t=0 Myr (the present day). 
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Extended Data Figure 3: PDF of the estimate of the number of supernovae required to power the Local 
Bubble’s expansion. The estimate is obtained by comparing the shell’s present-day momentum to the average 
momentum injected by supernovae. 
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Extended Data Figure 4: Analysis of the stellar tracebacks of the UCL and LCC clusters, whose progenitors 
were likely responsible for the supernovae that created the Local Bubble. The scatter points indicate the 
positions of the current cluster members of UCL and LCC, which are colored as a function of time (spanning the 
present day in pink to 30 Myr ago in black). Panel a: Using Hipparcos data and adopting a solar peculiar motion 
(Uo, Vo, We) = (10.0, 5.2, 7.2) km/s*, prior literature®’ find that UCL and LCC were born outside the Local Bubble 
(black trace’) 15 Myr ago and only entered its present-day boundary in the past 5 Myr (reproduced from Figure 
6, The search for the origin of the Local Bubble redivivus, Fuchs et al. 2006. Mon. Not. R. Astron. Soc. Volume 373, Issue 
3, pp. 993-1003). Panel b: We attempt to reproduce the results from prior literature®’ using the same data and 
solar motion, but are unable to account for the curvature of the traceback, finding the UCL and LCC formed just 
inside its northern boundary 15 Myr ago. Panel c: Using a different value for the solar motion, (Uo, Vo, We) = 
(10.0, 15.4, 7.8) km/s*! but the same Hipparcos data, we find that UCL and LCC were born near the center of the 
Local Bubble. Panel d: Finally, using updated Gaia data but the same adopted solar motion used in panel c. (Uo, 


Vo, Wo) = (10.0, 15.4, 7.8) km/s“, we also find that UCL and LCC were born near the center of the bubble, given 


an updated model for its surface.’ 


